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Abstract 

A  model  for  bottom  backscattering  strength  in  the  frequency  range  100  Hz 
to  10  kHz  is  presented.  This  model  generalizes  an  earlier  low-frequency 
model  intended  for  the  smaller  frequency  range  100-1000  Hz.  The  new  model 
matches  the  predictions  of  an  existing  high-frequency  model  (10-100  kHz) 
but  offers  new  insights  into  the  effect  of  fine-scale  layering.  Scattering  from 
both  sediment  volume  inhomogeneity  and  interface  roughness  is  included  in 
the  model.  For  silt  and  clay  sediments,  sediment  volume  scattering  is  usually 
the  dominant  process.  The  model  requires  profile  functions  for  average  sound 
speed,  mass  density,  sound  absorption,  and  volume  scattering  strength  in  the 
sediment.  In  addition,  two  parameters  describing  the  spectrum  of  interface 
roughness  are  required.  Use  of  the  small-roughness  perturbation  approxi¬ 
mation  prevents  application  of  the  model  to  extremely  rough  bottoms  (e.g., 
gravel  and  rock). 

Preliminary  comparisons  with  data  show  that  the  model  can  reproduce  an 
interesting  frequency  dependence  often  seen  in  bottom  backscattering.  This 
includes  a  decrease  in  scattering  strength  with  increasing  frequency  in  the  100- 
1000  Hz  range,  an  increase  as  frequency  approaches  10  kHz,  and  approximate 
constancy  or  a  slow  rise  above  10  kHz.  The  increase  in  sound  speed  with  depth 
typically  observed  in  sediments  causes  upward  refraction  which,  in  some  cases, 
leads  to  enhancements  of  scattering  by  6-8  dB. 
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1  INTRODUCTION 

This  report  describes  a  model  for  bottom  acoustic  backscattering  in  the  fre¬ 
quency  range  100  Hz  to  10  kHz.  This  model  is  a  generalization  of  a  model  covering 
the  low-frequency  range  100  Hz  to  1  kHz  [1,  2].  The  older,  low-frequency  model 
was  itself  an  adaptation  of  a  high-frequency  model  covering  the  frequency  range  10- 
100  kHz  [3].  The  present  work  bridges  the  1-10  kHz  gap  between  the  low-frequency 
and  high-frequency  regimes  in  a  manner  consistent  with  data  and  with  both  of  the 
older  models.  Rather  than  making  phenomenological  modifications  to  the  older 
model,  we  have  generalized  it  to  embody  both  large-  and  small-scale  environmental 
features.  This  more  general  model  contains  the  older  models  as  limiting  cases,  but 
also  allows  more  realistic  modeling  of  bottom  acoustic  properties. 

Like  the  earlier  models,  the  present  model  predicts  bottom  scattering  strength  in 
terms  of  sediment  acoustic  properties  (including  volume  scattering)  and  parameters 
for  interface  roughness.  Earlier  model/data  comparisons  have  shown  that  sediment 
volume  scattering  is  a  major  contributor  to  bottom  backscattering  [2,  4,  5],  and  it 
is  a  major  part  of  the  present  work.  In  the  high-frequency  model,  sediment  vol¬ 
ume  scattering  is  treated  in  an  approximation  in  which  gradients  in  average  sound 
speed,  absorption,  and  mass  density  are  neglected.  The  resulting  bottom  scattering 
strength  is  almost  independent  of  frequency,  provided  certain  assumptions  (to  be 
discussed  later)  are  made.  This  is  in  agreement  with  available  high-frequency  data 
for  mud  bottoms  [4,  5,  6]  where  sediment  volume  scattering  is  expected  to  dominate 
scattering  due  to  interface  roughness.  In  the  low-frequency  model,  average  sound 
speed  increases  with  depth  in  the  bottom  on  scales  of  10-100  m.  The  resulting 
upward  refraction  has  strong  effects  in  scattering  predictions.  The  most  impor¬ 
tant  refractive  effect  is  focusing  with  an  attendant  increase  in  scattering  strength. 
Figure  1  shows  this  increase  in  bottom  scattering  strength  for  frequencies  below 
1000  Hz,  an  effect  that  is  often  observed  in  real  data.  Wideband  data  often  exhibit 
a  U-shaped  curve  of  scattering  strength  vs  frequency,  with  scattering  strength  show¬ 
ing  a  minimum  somewhere  between  1  kHz  and  5  kHz.  This  minimum  is  above  the 
frequency  range  of  the  low-frequency  model,  but  if  the  model  is  used  in  this  range, 
it  produces  a  flat  curve  as  illustrated  in  Fig.  1. 

Although  the  high-frequency  model  ignored  gradients  in  sediment  properties, 
there  is  little  doubt  that  gradients  are  relevant  at  high  frequencies  as  well.  In  fact, 
the  high-frequency  model  uses  adjusted  values  of  sediment  sound  speed  and  density 
[3]  to  partially  account  for  the  steep  gradients  in  the  upper  few  centimeters  of  the 
bottom,  where  there  is  a  tendency  for  sediment  properties  to  approach  those  of 
water.  The  approach  we  take  in  the  present  work  is  to  allow  more  general  profiles 
of  sediment  acoustic  properties  with  gradients  on  both  small  scales  and  large  scales. 
In  this  approach,  the  transition  between  low-  and  high-frequency  models  occurs 
naturally. 
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Fig.  1.  Frequency  dependence  of  low-frequency  model  at  a  grazing  angle  of  45°  compared 
with  historical  data.  Three  values  of  the  initial  sound  speed  gradient,  g,  are  used.  As  the 
gradient  is  reduced,  the  model  curve  flattens  and  becomes  approximately  independent  of 
frequency. 


Figure  2  shows  typical  sound-speed  profiles  employed  in  the  low-  and  high- 
frequency  models,  along  with  a  more  realistic  profile.  The  term  realistic  should 
not  be  interpreted  to  mean  “typical,”  as  there  is  great  variety  in  the  geological  struc¬ 
ture  of  the  ocean  bottom.  The  more  realistic  profile  of  Fig.  2  includes  a  slow  upper 
layer  about  10  cm  thick.  The  upper  5-20  cm  of  the  bottom  is  often  a  region  of  strong 
mixing  due  to  bioturbation  (animal  activity)  and  hydrodynamic  activity,  which  can 
loosen  fine  sediments.  As  mentioned  earlier,  the  effects  of  this  layer  were  approx¬ 
imated  in  the  high-frequency  model  by  the  simple  expedient  of  reducing  sediment 
sound  speed  and  density  relative  to  conventional  values  obtained,  for  example,  from 
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Fig.  2.  Idealized  bottom  sound  speed  profiles  (dashed  curves)  employed  in  the  high- 
frequency  and  low-frequency  models  compared  with  a  more  realistic  profile  (solid  curve). 
The  sediment  begins  at  zero  depth.  The  water  sound  speed  (not  shown)  may  be  slightly 
larger  or  smaller  than  the  surficial  value  of  sediment  sound  speed. 


the  Hamilton-Bachman  regression  relations  [7].  The  adjusted  values  were  taken  to 
be  independent  of  depth  on  the  assumption  that  acoustic  penetration  would  not 
exceed  the  depth  of  the  upper  layer.  Deeper  in  the  sediment,  sound  speed  increases 
owing  to  consolidation.  In  the  older  low-frequency  model,  this  increase  is  gradual 
near  the  interface  (typically  about  1  m/s  per  meter  of  depth),  but  accelerates  and 
becomes  singular  at  finite  depth.  More  realistically,  the  sound  speed  at  depth  would 
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not  increase  past  the  value  expected  for  fully  consolidated  sediment  (e.g.,  mudstone) 
or  for  rock  (about  2000-6500  m/s)  [8].  Thus,  a  more  realistic  sound-speed  profile, 
as  indicated  in  Fig.  2,  includes  a  5-  to  20-cm-thick  slow  layer  near  the  interface,  and 
then  increases  with  depth  on  a  scale  of  order  100  m,  finally  reaching  a  maximum 
value.  In  the  region  where  the  maximum  sound  speed  is  attained,  the  bottom  ma¬ 
terial  is  fully  consolidated  sediment  or  rock  and  it  can  support  shear  waves.  Our 
present  model  still  does  not  treat  shear  effects,  which  are  most  likely  to  be  important 
only  when  the  sediment  cover  is  thin  or  the  frequency  is  very  low.  Unconsolidated 
sediments  do  support  shear  waves,  but  shear  effects  are  expected  to  be  very  small 
in  this  case  [9]. 

The  generalized  model  solves  a  problem  arising  from  the  artificial  nature  of  the 
sound  speed  profile  in  the  older  low-frequency  model.  In  that  model,  the  sound 
speed  profile  becomes  singular  (infinite)  at  a  finite  depth  in  the  sediment.  While 
this  singularity  is  usually  below  the  penetration  depth  of  acoustic  energy,  it  some¬ 
times  is  not,  particularly  for  steep  grazing  angles  and  at  the  low  frequencies  near 
100  Hz.  Mourad  and  Jackson  [2]  have  noted  that  this  singularity  causes  an  unusually 
strong  oscillation  in  the  bottom  scattering  strength  at  large  grazing  angles.  Such 
oscillations  are  an  expected  consequence  of  the  refracted  up-going  wave,  but  the 
singularity  gives  an  unrealistically  large  amplitude  for  the  up-going  wave  at  large 
grazing  angles.  This  large- amplitude  up-going  wave  also  causes  a  relatively  steep 
increase  in  scattering  strength  at  the  larger  grazing  angles  (roughly,  45°  to  85°).  It 
will  be  seen  that  this  steep  increase  is  absent  with  more  realistic  profiles. 

This  report  is  organized  as  follows.  Section  2  treats  propagation  of  sound  in 
the  bottom  and  volume  scattering  from  scatterers  or  inhomogeneities  embedded  in 
the  bottom  sediment.  It  is  mainly  this  portion  of  the  model  that  is  significantly 
improved  relative  to  earlier  versions,  but  there  is  also  a  related  improvement  in  the 
treatment  of  scattering  from  interface  roughness.  The  volume-scattering  component 
of  the  model  is  illustrated  in  Section  2  by  several  numerical  examples.  Section  3  gives 
a  brief  account  of  the  interface  roughness  scattering  portion  of  the  model,  which  has 
been  generalized  to  include  effects  of  gradients  in  the  sediment.  These  effects  are  also 
illustrated  by  numerical  examples.  Section  4  comments  upon  the  combined  model 
for  bottom  backscattering  and  discusses  dimensionless  sediment  profile  functions. 
Finally,  conclusions  and  recommendations  for  future  work  are  given  in  Section  5. 
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2  BOTTOM  PROPAGATION,  REFLECTION, 
AND  VOLUME  SCATTERING 

A  description  of  propagation  in  a  stratified  ocean  bottom  is  required  for  com¬ 
puting  both  sediment  volume  scattering  and  interface  roughness  scattering.  The 
earlier  low-frequency  model  [1,  2]  employed  a  linear  profile  for  the  squared  wave- 
number,  which  is  specific  to  Airy-function  solutions.  The  main  changes  in  the  present 
model  with  respect  to  this  earlier  work  are  as  follows: 

1.  The  wave  equation  is  now  solved  by  a  general  numerical  technique  for  smooth 
profiles  of  sediment  properties,  rather  than  for  the  specific  case  having  an  Airy- 
function  result.  This  allows  sediment  profiles  to  be  less  specific  than  before,  yet  still 
allows  checks  against  specific  prior  results. 

2.  Density-gradient  effects  are  treated. 

3.  Complex  reflection  coefficient  values  are  no  longer  inconsistent  with  the 
evaluation  of  sediment  propagation.  This  improves  upon  the  Rayleigh  impedance- 
mismatch  reflection  coefficient,  which  is  strictly  appropriate  only  for  a  semi-infinite 
sediment  of  depth-independent  properties. 

4.  Volume  scattering  and  absorption  within  the  ocean  bottom  are  depth  de¬ 
pendent. 

2.1  Bottom  Propagation 

Consider  horizontally  stratified  (fluid)  bottom  sediments  in  contact  with  isove¬ 
locity  seawater  at  a  horizontal  interface.  Shear  waves  are  not  treated  here.  The 
sediment  properties  of  sound  speed,  absorption,  and  density  p{z )  may  vary  with  the 
vertical  coordinate  z  in  the  sediments.  In  this  case  [10]  the  field  that  obeys  the  wave 
equation,  and  (for  a  given  acoustic  frequency)  also  the  Helmholtz  equation,  is  not 
the  acoustic  pressure  p,  but  rather  p/ ^ fp . 

For  incident  plane  waves  of  frequency  /  and  grazing  angle  6,  the  separation-of- 
variables  technique  yields 

P(x,z,t)lJ^7)  =  X(x)Z(z)e-2”'-,  (1) 

then  fix  the  value  of  a  0-dependent  separation  constant  k \  and  solve  the  remaining 
Helmholtz  equation  for  the  2-dependent  factor  Z(z).  That  is 

(2) 
(3) 


X"(x)  +  k2xX(x)  =  0  , 
X(x)  =  eik*x  , 
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and 

Z"(z)  +  Kl(z)Z(z)  =  0  .  (4) 

The  squared  wavenumber  field  needed  for  the  latter  equation  is  then 

Kl(z )  =  k\z)  -kl  +  i p"(z)/p(z)  -  \[p'(z)lp(z)}\  (5) 

with 

kx  =  2nf  cos  6  /  cw  .  (6) 

Here  the  ^-independent  constant  kx  is  the  (real)  horizontal  component  of  the  wave- , 
number,  conveniently  evaluated  here  at  z  —  0+  (the  sediment-water  interface  is  at 
^  =  0).  The  sound  speed  in  water  is  denoted  cw ,  and  k2(z)  is  the  (complex)  squared 
wavenumber  in  the  sediment.  The  imaginary  part  of  k(z)  allows  for  absorption  of 
sound  in  the  sediment.  The  density  terms  that  are  included  here  were  developed 
elsewhere  [10].  It  is  implicit  in  Eq.  (5)  that  profiles  of  sound  speed  and  absorption  are 
available  (to  construct  k2(z))  and  that  profiles  of  density  and  its  first  two  derivatives 
are  also  available. 

The  second-order  Helmholtz  equation  for  the  complex  function  Z(z)  is  obtained 
numerically  by  solving  the  equivalent  set  of  four  coupled  first-order  ordinary  differ¬ 
ential  equations  obeyed  by  functions  that  represent  the  separated  real  and  imaginary 
parts  of  Z(z)  and  Z'(z )  [11].  After  normalizing  Z(z)  for  a  plane  incident  sound  wave 
for  which  p/^/p  has  unit  amplitude,  the  bottom  backscattering  differential  cross  sec¬ 
tion  per  unit  area  and  per  unit  solid  angle  is  given  by  the  following  integration  over 
the  volume  scatterers  [1,  2]: 

roo 

crbv(6)  =  I  dz  av(z )  \Z(z)j4.  (7) 

Jo 

The  bottom  scattering  strength  as  conventionally  defined  is 

BSSV  =  10  log10  <rbv{0),  (8) 

where  the  V  indicates  that  we  are  considering  only  the  contribution  of  sediment 
volume  scattering  to  the  bottom  scattering  strength.  In  Eq.  (7),  av(z)  is  the  volume- 
specific  differential  scattering  cross  section  for  the  sediments,  i.e.,  10  log10  &v(z)  is 
the  “volume  scattering  strength”  [12]. 

2.2  Numerical  Solution 

The  numerical  solution  proceeds  upward  to  the  water /sediment  interface  from 
a  chosen  starting  depth  in  the  sediment.  Initial  values  for  Z  and  Z'  at  the  starting 
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depth  are  chosen  consistent  with  the  tunneling  or  oscillatory  character  of  Z(z)  there, 
that  is,  the  negative  or  positive  sign  of  the  real  part  of  K* ,  respectively. 

For  the  usual  case  where  relatively  small  grazing  angles  are  of  interest,  an  in¬ 
creasing  sediment  sound  speed  produces  a  turning  point  (or  vertex)  in  the  raytrac¬ 
ing  limit.  The  corresponding  wave  solution  is  exponentially  small  below  the  turning 
point  depth.  In  wave  theory,  this  is  the  depth  below  which  the  real  part  of  is 
negative.  An  unnormalized  wave  solution  having  arbitrary  amplitude  is  started  from 
well  into  this  evanescent  region  as  for  a  decaying  exponential,  so  that  derivative  and 
value  are  connected  according  to 

Z'/Z  =  \K,\  .  (9) 

The  amplitude  of  this  wave  is  fixed  by  a  normalization  procedure  defined  below. 

In  steep- angle  cases,  rays  need  never  vertex,  but  can  instead  continue  to  arbi¬ 
trarily  great  depths  and  eventual  absorption.  The  wave  solution  is  then  an  entirely 
downward- directed  field  at  depths  below  all  changes  in  sediment  properties.  Then 
the  solution  is  started  as  for  a  traveling  wave,  so  that 

Z'/Z  =  -i  I Kz\  .  (10) 

In  all  cases  the  integration  proceeds  upward  from  the  starting  depth  until  the 
water /sediment  interface  is  reached.  There  the  solution  is  normalized  by  forcing 
continuity  across  the  interface  for  both  the  acoustic  pressure  p(z)  =  Z(z)  \J p{z) 
and  the  ^-derivative  quantity  p\z)/p(z). 

The  complex  reflection  coefficient  for  the  water /sediment  interface  is  the  ratio  of 
upward  and  downward  wave  amplitudes  in  the  wave  solution  on  the  water  side  of  the 
interface.  In  this  wave-theory  treatment,  the  reflection  coefficient  depends  strongly 
on  angle  and  frequency,  unlike  the  simpler  Rayleigh  impedance-mismatch  reflection 
coefficient  appropriate  to  a  uniform  sediment  having  depth-independent  sound  speed 
and  density.  These  improved  values  for  the  reflection  coefficient  apply  even  in  the 
Airy-function  (it2- linear-profile)  limit.  This  improvement  influences  estimates  of 
scattering  by  interface  roughness,  which  is  sensitive  to  the  relative  amplitudes  of 
upward-directed  and  downward- directed  fields  at  the  interface-roughness  scatterers. 
As  will  be  seen  in  the  following  section,  the  influence  on  interface  scattering  can  be 
very  strong. 

The  normalized  wave  solution  is  obtained  after  jumping  the  wave  solution  into 
the  water  across  the  water /sediment  interface  and  renormalizing  for  unit  amplitude 
of  the  downward-directed  incident  plane- wave  component  of  the  wave  solution  above 
the  sediment.  In  the  water  above  2  =  0,  let 

Z(z)  =  De+ik*z  +  Ue~ik’z  ,  (11) 
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with 

kz  =  — 2tt/  sin  0  /  cw  . 

(12) 

Approaching  z  =  0  from  the  water  side,  we  then  have 

(Z  /Z)0+  -  ikz  +  . 

(13) 

The  complex  reflection  coefficient  is  just 

(14) 

In  the  evaluation  of  roughness  scattering,  the  above  result  also  simplifies  the  fol¬ 
lowing  expression  that  appears  in  the  roughness-scattering  portion  of  the  model.  We 
note  that  the  normalized  admittance  presented  by  the  interface  to  a  plane  acoustic 


wave  is 


1  -  R(0)  (Z'/Z)0+ 


(15) 


1  +  R(6)  ikz 

The  unnormalized  solution  is  just  divided  everywhere  by  the  complex  amplitude  D 
of  its  downgoing  component  in  the  water  to  produce  the  normalized  Z(z)  solution 
for  unit  incident-wave  amplitude.  This  normalization  step  is  not  needed  for  the 
complex  reflection  coefficient  U/ D,  but  it  is  needed  for  sediment  scattering-strength 
integrals. 

Bottom  reflection  loss  is  given  by 


BRL  =  -20  log10  \U/D\ 


(16) 


This  is  a  positive  decibel  equivalent  of  the  complex  reflection  coefficient  and  is  shown 
along  with  the  scattering  strength  BSSV  in  some  of  the  examples  that  follow. 

The  following  examples  illustrate  some  of  the  features  of  the  generalized  treat¬ 
ment  of  sediment  volume  scattering.  Low-frequency  examples  are  considered  for 
the  most  part,  as  propagation  effects  are  most  obvious  below  about  1  kHz.  It  is 
important  to  note,  however,  that  the  generality  of  the  model  makes  it  possible  to 
treat  the  effects  of  sediment  layering  at  both  large  and  small  scales,  and  hence  at 
both  low  and  high  frequencies.  Figure  3  gives  the  sound  speed  profiles  used  in  the 
following  examples. 
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Sound  Speed  (m/s) 


Fig.  3.  Sediment  sound  speed  profiles  used  in  illustrative  examples. 


2.3  A;2-Linear  Profile 


Figures  4  and  5  show  calculations  of  BSSV  (from  sediment  volume  scattering 
only — no  roughness)  that  check  known  cases  from  prior  work,  except  that  the  wave 
solution  is  now  obtained  by  a  more  general  numerical  technique.  In  this  prior  work, 
k2(z )  was  taken  to  have  a  linear  dependence  on  depth  [1,  2],  for  which  the  field  Z(z) 
is  given  by  an  Airy  function.  The  parameters  and  profile  functions  appropriate  to  a 
comparison  with  Mourad  and  Jackson  are 


Frequency,  / 
Water  sound  speed,  cw 
Surficial  sediment  sound  speed,  c\ 
Initial  slope  of  sound  speed  profile,  \dc/dz\0 

Loss  parameter,  S 


500  Hz  and  100  Hz 
1530.0  m/s 

1510.11  m/s  =  cw  x  0.987 

1.0  s'1 

0.00164 
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Squared  wavenumber,  k2(z)  — 
Real  part  of  surficial  wavenumber,  = 
Depth  of  singularity,  2wan  = 
Mass  density,  p  = 
Absorption  coefficient,  ai  = 
Volume  scattering  parameter,  a2  = 
Sediment  volume  scattering  strength,  av  = 


k\  (1  +  f£)2  (1  -  z/zy,^) 
2?r  f/c! 

—755.1  m 

1.72  x  1026  kg/m3 

407 v8f  /  (  Ci  ln(10)  ) 

0.0002 

ai  a2  • 


The  parameter  2waU  is  the  depth  at  which  the  sound  speed  becomes  infinite 
for  this  special  profile.  It  is  not  an  independent  parameter  and  can  be  found  in 
terms  of  c\  and  \dc/dz\o.  The  parameter  <y\  is  the  sound  absorption  coefficient  (in 
decibels/unit  length)  just  below  the  sediment /water  interface,  and  <r2  [1,  2,  3]  is  a 
convenient  dimensionless  form  for  the  volume  scattering  strength,  as  discussed  in 
Section  4.  The  parameters  and  results  here  are  the  same  as  those  in  prior  calculations 


curve  is  computed  for  the  gradient-free  case. 
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GRAZING  ANGLE  (deg) 

Fig.  5.  As  in  Fig.  4  except  the  frequency  is  100  Hz. 


for  this  Airy-function  case  (see  Fig.  3  of  Ref.  2).  This  is  a  significant  check  of  the 
present  numerical  methods.  The  500-Hz  BSSV  values  show  little  of  the  oscillatory 
interference  structure  vs  grazing  angle  that  is  so  apparent  at  100  Hz  (see  Fig.  4). 
This  effect  comes  from  the  higher  sound  absorption  at  the  higher  frequency  and 
the  consequently  weaker  upward-directed  wave  reflecting  off  the  deep  barrier  in  the 
sediment.  The  absorption  of  sound  at  100  Hz  is  so  low  that  interference  of  the 
upward  and  downward  waves  in  the  sediment  is  important  [2]. 

Also  shown  with  Figs.  4  and  5  are  the  corresponding  values  for  bottom  reflection 
loss,  BRL:  and  the  smooth,  frequency-independent  Rayleigh  values  calculated  from 
the  impedance  mismatch  at  the  interface  between  sediment  and  water.  The  BRL 
curves  emphasize  the  phenomena  that  have  arisen  at  low  frequencies  as  a  result  of 
the  combination  of  upward  refraction  and  small  absorption.  The  500-Hz  results  of 
Fig.  4  show  that  for  most  angles  BRL  can  be  approximated  by  the  Rayleigh  reflection 
coefficient  computed  from  the  impedance  mismatch.  The  deepest  sediments  are  not 


TM1-94  11 


UNIVERSITY  OF  WASHINGTON  -  APPLIED  PHYSICS  LABORATORY, 


important  when  sound  absorption  is  high.  The  100-Hz  results  of  Fig.  5,  however, 
show  the  systematically  lower  reflection  losses  that  correspond  to  low-frequency 
sound  reaching  and  returning  from  deep  sediments.  Comparing  the  peaks  in  BSSV 
and  BRL  shows  an  obvious  correlation  between  high  BSSV  and  high  BRL.  That 
is,  greater  sediment  backscattering  goes  hand  in  hand  with  greater  penetration  of 
sound  into  the  sediments  below  the  interface.  Deep  sediments  are  important,  and 
different  possible  treatments  of  them  are  distinguishable  at  low  frequencies. 

2.4  fc2-Linear  Profile  with  Maximum  Sound  Speed 

Figure  6  shows  the  effect  of  clipping  the  sound  speed  at  a  maximum  value 
rather  than  letting  it  become  infinite  at  z  —  2wau  and  then  imaginary  at  greater 
depths,  as  it  does  for  the  Airy-function  limit.  The  parameters  for  the  calculations 
in  Fig.  6  are  the  same  as  those  for  the  100-Hz  calculations  of  Fig.  5,  except  that  for 


Fig.  6.  As  in  Fig.  5  except  the  sound  speed  was  held  constant  at  depths  greater  than 
the  depth  at  which  it  is  twice  the  surficial  value. 
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these  calculations  the  sediment  sound  speed  was  only  allowed  to  double  in  value  with 
increasing  depth.  At  greater  depths,  it  was  then  held  constant.  Near  90°  grazing 
angle,  these  low-frequency  results  are  quite  different  from  the  Airy  case,  as  expected 
from  the  fact  that  rays  would  no  longer  vertex  at  an  unphysical  barrier  but  would 
instead  continue  deeper  and  deeper  until  absorbed.  With  the  absence  of  an  upward 
wave  reflected  off  the  deep  barrier  at  such  steep  incidence,  the  interference  character 
in  BSSV  has  disappeared.  The  impedance-mismatch  values  for  BRL  have  also 
become  a  good  approximation.  At  shallower  angles,  however,  the  unclipped  profile 
does  very  well,  giving  essentially  identical  results  (compare  Figs.  5  and  6).  This  also 
accords  with  ray  theory,  since  rays  can  only  reach  the  changed  part  of  the  sediment 
profile  for  the  steepest  grazing  angles.  The  most  striking  feature  of  Fig.  6  is  the 
approximately  6-dB  drop  in  BSSV  for  the  steep  grazing  angles  where  there  is  no 
upward-refracted  wave.  While  this  is  not  a  very  realistic  situation,  it  illustrates 
the  strong  scattering  enhancements  that  occur  as  a  result  of  upward  waves.  From  a 
ray  viewpoint,  this  enhancement  corresponds  to  superposition  of  the  four  round-trip 
multipath  combinations  provided  by  upward  refraction.  This  can  give  a  6-dB  effect, 
although  the  upward- refracted  wave  increases  average  sound  intensity  at  the  volume 
scatterers  by  3  dB  at  most  (when  the  amplitude  of  the  upward  wave  is  comparable 
to  that  of  the  downward  wave). 

At  higher  frequencies,  even  near-vertical  backscattering  is  not  sensitive  to  the 
presence  or  absence  of  a  barrier  at  z  =  2wau  (e.g.,  see  Fig.  4).  The  higher  sound 
absorption  makes  deep-sediment  details  irrelevant.  The  sediments  near  the  water 
interface  then  dominate. 


2.5  Linear  Sound  Speed  Profile 


Figure  7  shows  100-Hz  BSSV  calculations  for  a  linear  sound-speed  profile. 
Here,  the  slope  of  the  linear  profile  is  the  same  as  the  initial  slope  in  the  linear-  k2 
Airy  case  of  Fig.  5,  which  corresponded  to  infinite  sound  speed  755.1  m  into  the 
sediment.  An  important  difference  between  these  cases  is  that  different  amounts  of 
the  sediment  are  strongly  ensonified.  In  particular,  the  linear-c(z)  case  does  not  have 
an  unphysical  barrier  for  vertical  ensonification.  Parameters  and  profile  functions 
used  for  Fig.  7  were 


Frequency,  / 
Water  sound  speed,  cw 
Surficial  sediment  sound  speed,  Cj 
Slope  of  sound  speed  profile,  \dc/dz\ 
Loss  parameter,  8 
Squared  wavenumber,  k2(z) 
Real  part  of  surficial  wavenumber,  k\ 


100  Hz 
1530.0  m/s 

1510.11  m/s  =  cw  x  0.987 
1.0  s"1 


0.00164 

k\  (1  +  i8)2  /  (1  +  z/zTet)2 
2tt//ci 
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Fig.  7.  As  in  Fig.  5  except  the  sound  speed  (rather  than  squared  wavenumber)  varies 
linearly  with  depth. 


Reference  depth,  zre f  =  —1510. 2  m 

Mass  density,  p  =  1.72  x  1026  kg/m3 
Volume  scattering  parameter,  ct2  =  0.0002  . 

The  volume  scattering  strength  is  derived  from  the  dimensionless  parameter  a2 
as  in  the  previous  examples.  The  parameter  zTe f  determines  the  slope  of  the  sound 
speed  profile.  The  low-angle  BSSV  results  for  the  linear-c(z)  case  of  Fig.  7  are 
similar  to  the  low-angle  results  in  the  Airy  case  of  Fig.  5.  However,  the  steeper  angles 
show  a  systematic  difference,  as  might  be  expected  from  a  raytracing  viewpoint. 
Once  again,  the  deep  details  of  the  sediments  have  become  significant  at  this  low 
frequency.  Absence  of  the  “wall”  has  reduced  the  magnitude  of  the  oscillations  in 
BSSV  as  a  function  of  grazing  angle  and  has  reduced  the  scattering  strength  at 
angles  greater  than  about  40°. 
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2.6  Linear  Sound  Speed  Profile  with  Maximum  Sound  Speed 

Figure  8  shows  results  calculated  for  yet  another  shape  of  the  c(z )  sound- 
speed  profile.  The  sound  speed  profile,  c(z),  is  again  linear  but  is  now  held  constant 
starting  at  the  depth  of  a  basement  under  the  sediment.  This  case  was  motivated  by 
a  published  sediment  raytracing  example  for  a  site  in  the  Caribbean  Sea  (see  Fig.  8 
of  Ref.  13).  Absorption  and  maximum  sound  speed  were  treated  as  in  the  prior 
examples  above,  but  with  the  altered  numerical  values  shown  below.  Otherwise, 
values  were  chosen  to  be  similar  to  the  raytracing  example  cited. 

Frequency,  /  =  100  Hz 
Water  sound  speed,  cw  =  1547.79  m/s 
Surficial  sediment  sound  speed,  C\  =  1508.17  m/s  =  cw  x  0.9744 
Slope  of  sound  speed  profile,  \dc/dz\  =  1.3  s-1 
Loss  parameter,  6  =  0.0015 


0  20  40  60  80 

GRAZING  ANGLE  (deg) 

Fig.  8.  Model  curves  for  reflection  loss  and  backscattering  strength  for  a  linear  sound- 
speed  profile  with  a  constant-speed  basement.  The  frequency  is  100  Hz. 
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Squared  wavenumber,  k2(z) 

Real  part  of  surficial  wavenumber,  ki 
Reference  depth,  zre { 
Basement  depth, 

Mass  density,  p 
Volume  scattering  parameter,  <r2 


k\  (1  +  iS)2  /  (1  +  zjzTe{)2 
(for  2  above  2base) 

2  JT//C! 

-1160.1  m 

—1097.2  m  =  —600  fathoms 
1.5  x  1026  kg/m3 
0.001. 


The  volume  scattering  strength  is  derived  from  the  dimensionless  parameter  <r2 
as  in  the  previous  examples. 


2.7  Frequency  Dependence  Due  to  Strongly  Scattering 
Upper  Sediment 

Figure  9  shows  total  BSS  values  for  several  provinces  from  the  Geddes  compi¬ 
lation  of  historical  data  [14].  The  figure  gives  backscattering  strength  vs  frequency 
for  a  grazing  angle  of  45°.  These  experimental  results  for  backscattering  have  a  U- 
shaped  frequency  dependence,  with  the  minimum  backscattering  in  our  frequency 
range.  This  U  shape  is  seen  in  the  data  from  several  but  not  all  of  the  provinces  in 
Ref.  14. 

To  generate  such  an  effect  in  the  simplest  possible  way,  it  was  found  sufficient  to 
increase  the  volume  scattering  parameter  cr2  for  a  thin  upper  portion  of  the  seabed 
sediments,  without  giving  its  value  any  frequency  dependence  or  a  more  detailed 
profile  function.  The  model  then  yields  a  fall  and  a  rise  in  backscattering  arising 
naturally  from  the  frequency  dependence  in  the  absorption  loss  (constant  6  vs  fre¬ 
quency).  Two  competing  effects  account  for  this.  First,  absorption  rises  linearly 
with  frequency  for  a  frequency-independent  6  value,  so  that  BSSV  values  initially 
fall  with  increasing  frequency  owing  to  supression  of  the  upgoing  wave.  (The  6 
constant  is  a  conversion  factor  giving  the  imaginary  part  of  the  wavenumber  from 
the  real  part,  which  rises  with  frequency,  so  that  absorption  also  rises.)  The  thin 
upper  sediment  layer  is  almost  irrelevant  at  low  frequencies  because  the  depth  of 
ensonification  is  much  greater  than  the  layer  thickness.  As  the  frequency  increases, 
the  upper  sediment  layer  becomes  important  compared  with  the  deeper  sediment 
as  increased  absorption  reduces  the  depth  of  penetration,  so  that  upper-sediment 
scattering  controls  the  BSSV  values  at  high  frequencies.  Finally,  at  frequencies 
such  that  the  penetration  depth  is  smaller  than  the  layer  thickness,  BSSV  becomes 
frequency  independent  as  the  linear  increase  in  volume  scattering  strength  with  fre¬ 
quency  (dictated  by  the  relation  given  earlier  between  the  volume  scattering  strength 
<jv  and  the  dimensionless,  frequency-independent  parameter  a2)  offsets  the  inverse 
frequency  dependence  of  the  penetration  depth. 
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200  500  1000  2000  5000  10000  20000  200  500  1000  2000  5000  10000  20000 

f  (Hz)  f  (Hz) 

Pr4  Lower  Continental  Rise  Pr7  Bermuda  Pedestal 


200  500  1000  2000  5000  10000  20000  200  500  1000  2000  5000  10000  20000 

f  (Hz)  f  (Hz) 


200  500  1000  2000  5000  10000  20000  200  500  1000  2000  5000  10000  20000 

f  (Hz)  f  (Hz) 


Fig.  9.  Compilation  of  historical  bottom  backscattering  data  by  Geddes  showing  fre¬ 
quency  dependence  at  a  grazing  angle  of  45°.  The  straight  lines  are  regression  fits. 


Figure  10  compares  calculations  without  and  with  such  an  increase  in  the  volume 
scattering  function  crv(z )  for  the  upper  layers,  and  the  resulting  U-shaped  frequency 
dependence  in  the  latter  case.  For  these  calculations,  we  have  chosen  the  same  k2- 
linear  profile  as  for  Figs.  4  and  5  but  have  varied  the  frequency  from  100-20,000  Hz. 
Here  the  conditions  are  the  same  for  both  the  lower  and  upper  curves,  except  that 
the  upper  curve  has  a  30-cm  thick,  15-dB  stronger  upper  layer. 
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Fig.  10.  Model  curves  showing  frequency  dependence  of  reflection  loss  and  backscatter- 
ing  with  an  upper  sediment  layer  exhibiting  strong  volume  scattering. 


A  comparison  between  the  layered  model  and  historical  data  at  grazing  angles 
near  45°  is  made  in  Fig.  11.  The  data  come  from  the  Lower  Continental  Rise,  which 
corresponds  to  Geddes’  Province  4.  Both  Geddes’  average  of  historical  data  and 
the  DREA  data  [14]  (which  are  part  of  the  average  shown  in  Fig.  10)  are  shown 
in  the  figure.  In  the  model  computation,  we  assumed  that  the  P-linear  profile  was 
applicable  with  the  following  parameters: 


Frequency,  / 
Water  sound  speed,  cw 
Surficial  sediment  sound  speed,  ci 
Initial  slope  of  sound  speed  profile,  |c?c/dz|o 

Loss  parameter,  8 
Squared  wavenumber,  k2(z) 
Real  part  of  surficial  wavenumber,  ki 
Depth  of  singularity,  zwau 
Mass  density,  p 
Volume  scattering  parameter,  a2 


=  500  Hz  and  100  Hz 
=  1530.0  m/s 

=  1510.11  m/s  =  cwx  0.987 
=  1.0  s"1 
=  0.00149 

=  k\  (1  +  iS )2  (1  -  zjzy, all) 

=  2tt//ci 
=  -755.1  m 
=  1.72  x  1026  kg/m3 
=  0.0003353  . 
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Fig.  11.  Comparison  of  frequency  dependence  of  data  from  the  Lower  Continental  Rise 
of  the  Northeast  Atlantic  with  model  curves. 


The  model  input  parameters  were  in  part  determined  from  core  samples  as  explained 
in  Ref.  [15].  Three  different  choices  were  made  for  the  strongly  scattering  upper 
layer:  (1)  a  thickness  of  1  m  and  a  volume  scattering  strength  that  is  elevated 
by  10  dB  relative  to  the  value  specified  above,  (2)  a  thickness  of  0.5  m  and  a 
volume  scattering  strength  elevated  by  10  dB,  and  (3)  a  thickness  of  0.3  m  and  a 
volume  scattering  strength  elevated  by  15  dB.  The  latter  two  choices  give  reasonable 
fits  to  the  data,  suggesting  a  layer  of  strong  scattering  that  is  about  30  cm  thick. 
Comparison  of  Figs.  10  and  11  shows  that  the  model  curve  of  Fig.  11  is  smoother  at 
high  frequencies.  This  is  because  we  employed  the  Airy  function  software  as  in  Ref. 
[1]  rather  than  the  numerical  propagation  code.  The  propagation  code  was  having 
numerical  difficulty  at  the  highest  frequencies,  but  this  was  evidently  not  due  to  any 
fundamental  limitation  of  the  code,  but  rather  to  the  fact  that  the  starting  point  of 
the  integration  was  taken  to  be  too  deep  in  comparison  to  the  depth  of  penetration. 
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3  ROUGHNESS-SCATTERING  COMPONENT 

The  following  development  is  as  previously  given  by  Mourad  and  Jackson  [3], 
except  for  three  changes.  The  most  significant  change  is  a  generalization  to  include 
the  above  effects  of  gradients  in  sediment  properties  in  the  roughness- scattering 
component  of  the  model.  In  addition,  an  algebraic  fit  to  the  Kirchhoff  integral  is 
abandoned  in  favor  of  numerical  evaluation  of  the  integral.  Finally,  the  composite¬ 
roughness  approximation  is  abandoned  in  favor  of  the  simpler  small-roughness  per¬ 
turbation  approach.  In  the  first  place,  we  have  found  that  the  particular  slope¬ 
averaging  method  used  in  our  earlier  models  could  not  deal  with  the  oscillations 
in  scattering  strength  seen  in  the  present  model.  Moreover,  the  difference  between 
the  composite-roughness  and  the  small-roughness  perturbation  approximations  was 
found  to  be  very  small  at  frequencies  below  10  kHz.  Given  the  same  bottom  pa¬ 
rameters,  the  present  model  will  substantially  match  the  the  predictions  of  the  older 
high-frequency  model  except  for  very  rough,  hard  bottoms  (gravel  and  rock),  where 
the  high-frequency  model  uses  an  empirical  expression  [3].  The  present  model  is 
inapplicable  in  these  situations. 

As  in  the  prior  high-frequency  model  [3],  bottom  relief  is  assumed  to  be  a 
Gaussian  random  process  having  the  following  two-dimensional  spectrum: 

IF(K)  =  w2/(I<h0y  .  (17) 

In  this  expression,  w2  and  7  are  model  parameters  for  the  strength  and  exponent  of 
the  roughness  spectrum,  and  h0  —  1  m  is  simply  a  reference  length.  Since  the  rough¬ 
ness  spectrum  depends  only  on  the  magnitude  of  the  two-dimensional  wavevector 
K  =  |K|,  it  can  be  seen  that  bottom  roughness  is  assumed  to  be  isotropic.  The  spec¬ 
trum  is  defined  for  both  positive  and  negative  arguments  and  is  normalized  so  that 
the  integral  over  any  finite  region  of  K-space  is  equal  to  the  mean-square  roughness 
due  to  those  Fourier  components  included  in  the  integral.  Random  processes  having 
power-law  spectra  such  as  Eq.  (17)  are  fractal  without  a  well-defined  correlation 
scale  or  rms  roughness.  Even  so,  finite  and  useful  results  can  be  obtained  provided 
the  exponent  satisfies  2  <  7  <  4. 

The  roughness  statistics  can  be  alternatively  described  by  the  “structure  func¬ 
tion”  [16],  itself  a  power  law: 

<  [/(Ro  +  R)  -  /(Ro)]2  >=  C2hR2a  ,  (18) 

where  /( R)  is  the  ^-coordinate  of  the  interface  at  horizontal  position  R  =  (x,y). 
The  function  /( R)  is  a  zero-mean  Gaussian  random  process.  The  structure-function 
parameters  are  related  to  the  spectrum  parameters  as  follows: 

o  =  |-l,  (19) 
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and 

2  2tt  w2  T{2  -  a)  2~2a 
h~  K  a(l  —  a)  r(l  +  a)  ' 

The  roughness  scattering  component  of  the  backscattering  cross  section  will  be 
written  as  the  sum  of  two  terms: 

<?br(0)  =  akr  (%(z)  +  <?pr  (#)[l  ~  ^(z)]  •  (21) 


The  first  term  in  this  equation  is  computed  using  the  Kirchhoff  approximation,  and 
the  second  is  computed  using  the  small-roughness  perturbation  approximation.  The 
function,  g(x),  interpolates  between  these  two  approximations. 

g(x)  =  •  (22) 

1  +  ex 

The  argument  of  g{x ), 


I  =  80  |  0.03  ]  . 


(23) 


is  chosen  such  that  the  Kirchhoff  approximation  is  used  for  angles  near  vertical  in¬ 
cidence  and  the  small-roughness  approximation  is  used  for  all  other  angles.  Specif¬ 
ically,  the  Kirchhoff  approximation  is  used  from  90°,  which  is  an  extremum,  down 
to  the  angle  at  which  the  Kirchhoff  cross  section  has  dropped  by  a  factor  of  0.03 
(about  —15  dB)  relative  to  this  extremum. 


As  in  Ref.  4,  the  Kirchhoff  cross  section  will  be  expressed  in  the  following  form: 

|R(t/2)|2 


<Tkr(0)  = 


87r  sin2  6  cos2  6 


f 

Jo 


0-qu- 


J0{u)  u  du 


(24) 


where 

q  =  Cl  sin2  6  cos-2“  6  2x~2a  .  (25) 

In  Eqs.  (24)  and  (25),  kw  is  the  (real)  wavenumber  in  water,  R(tt/2)  is  the  complex 
plane- wave  reflection  coefficient  for  a  flat  interface  separating  water  and  sediment, 
here  evaluated  at  normal  incidence,  and  J0(u )  is  the  zeroth-order  Bessel  function  of 
the  first  kind. 


A  significant  departure  is  made  from  previous  versions  of  the  APL-UW  bottom 
backscattering  model,  in  which  gradients  were  ignored  and  the  Rayleigh  (or  Fresnel) 
reflection  coefficient  was  used.  To  account  for  gradients  in  the  sediment  properties, 
the  reflection  coefficient,  R{6 ),  is  now  computed  numerically,  as  developed  in  Sec¬ 
tion  2.2  above,  using  the  profiles  for  sound  speed,  density,  and  absorption  coefficient. 
Equation  (24)  is  also  evaluated  numerically,  as  in  Ref.  4.  This  improves  upon  the  al¬ 
gebraic  fit  to  Eq.  (24)  used  by  Mourad  and  Jackson  [3].  In  this  numerical  evaluation, 
we  use  an  algorithm  provided  by  C.  deMoustier  (private  communication). 
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In  the  small-roughness  approximation,  first-order  perturbation  theory  is  used 
to  compute  scattering  from  the  small-scale  roughness  of  the  interface.  Following 
Winebrenner’s  treatment  of  the  electromagnetic  scattering  problem,  Moe  and  Jack- 
son  [17]  have  obtained  the  following  expression  for  the  acoustic  backseat tering  cross 
section  of  a  randomly  rough  two-fluid  boundary: 


<v(*)  = 


^\1  +  R{6)\4  W(2kwcos6) 

x  |1  —  —  +  (1  —  — )  cos2  0  +  (pi  -  1) 
Pi  Pi 


i-R(d)i2 

i  +  R{0)- 


sin20  |2  .  (26) 


In  this  equation,  pi  is  the  ratio  of  sediment  density  immediately  below  the  interface 
to  water  density.  Similarly, 

k  =  ki  (1  -f  i8)/kw  (27) 


is  the  ratio  of  the  complex  wavenumber  in  the  sediment  immediately  below  the  in¬ 
terface  to  the  real  wavenumber  in  the  water.  Expression  (26)  is  a  rearrangement  of 
the  conventional  perturbation-theory  result  [3,  4,  18].  In  this  rearrangement  the  re¬ 
flection  coefficient  R(0)  for  the  (flat)  water/sediment  interface  appears  prominently. 
The  conventional  result  is  recovered  if  the  Rayleigh  reflection  coefficient  is  used,  and 
the  more  general  result  obtains  if  the  reflection  coefficient  is  computed  numerically, 
including  gradients. 

Figures  12a  and  12b  illustrate  the  effect  of  sound  speed  gradients  on  roughness 
scattering.  The  P-linear  profile  is  used  again,  allowing  computation  of  the  reflection 
coefficient  in  terms  of  Airy  functions.  The  relevant  parameters  are 


Frequency,  /  = 
Water  sound  speed,  cw  — 
Surficial  sediment  sound  speed,  C\  = 
Slope  of  sound  speed  profile,  \dc/dz\0  = 
Loss  parameter,  8  = 
Squared  wavenumber,  k2(z)  = 
Real  part  of  surficial  wavenumber,  kx  = 
Depth  of  singularity,  zwaa  = 
Mass  density,  p  = 
Volume  scattering  parameter,  cr2  = 
Spectral  exponent,  7  = 
Spectral  strength,  ic2  — 


300  Hz,  1000  Hz 
1545.0  m/s 

1483.2  m/s  =  Cw  x  0.96 

1.3  s"1 

0.0001445 

k\  (1  -+•  i6 )2  (1  -  z  j 2wau) 
2?r  f  I  C\ 

— 570.5  m 

1.548  x  1026  kg/m3 
0.0 
3.25 

0.0001636  m4. 


The  features  of  the  roughness-scattering  component  of  the  scattering  strength 
are  similar  to  those  seen  for  the  volume-scattering  component.  The  scattering 
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strength  exhibits  strong  oscillations,  and  upward  refraction  causes  an  enhancement 
of  the  backscattering  strength.  As  in  the  volume-scattering  cases  seen  earlier,  the 
oscillations  and  enhancement  weaken  as  frequency  increases,  owing  to  absorption 
which  suppresses  the  upgoing  wave.  The  oscillations  show  slow  variations  in  strength 
as  a  function  of  grazing  angle.  This  is  due  to  use  of  composite  roughness  slope  av¬ 
eraging  as  defined  by  Mourad  and  Jackson  [3].  Thus  the  computations  leading  to 
Fig.  12  depart  slightly  from  the  model  defined  in  this  section,  in  which  composite 
roughness  averaging  is  omitted.  The  sound-speed  profile  used  here  does  not  con¬ 
tain  the  type  of  small-scale  structure  that  would  be  important  at  high  frequencies. 
Bioturbation  causes  strong  gradients  in  sediment  acoustic  properties  on  scales  of 
5-20  cm  [5]  which  could  cause  upward  refraction  at  high  frequencies.  The  present 
model  provides  a  means  for  investigating  such  effects. 


Fig.  12a.  Model  curves  for  backscattering  strength  due  to  interface  roughness  with  a 
gradient  (solid  curve)  and  without  a  gradient  (dotted  curve),  300  Hz. 
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Fig.  12b.  Model  curves  for  backscattering  strength  due  to  interface  roughness  with  a 
gradient  (solid  curve)  and  without  a  gradient  (dotted  curve),  1000  Hz. 
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4  COMMENTS 

Sections  2  and  3  defined  the  methods  for  computing  the  backscattering  cross 
section  for  scattering  from  sediment  volume  inhomogeneities,  and  <7^(0),  for 

scattering  from  the  rough  interface.  Bottom  backscattering  strength  BSS  is  defined 
as  the  backscattering  cross  section  per  unit  area  per  unit  solid  angle,  expressed  in 
decibels.  This  is  equivalent  to  the  standard  definition  given  by  Urick  [12].  Combining 
the  two  components,  we  will  write  the  total  backscattering  strength  in  the  form 

BSS  =  Sb(0)  =  10  log  jo  K(0)  +  crbr(0)\  .  (28) 

The  bottom  material  is  treated  as  a  fluid  supporting  compressional  waves  only. 
Consequently,  the  sediment  properties  of  interest  correspond  to  those  of  the  APL- 
UW  high-frequency  model  [3],  except  that  the  sediment  properties  density,  sound- 
speed,  absorption,  and  volume  scattering  strength  are  allowed  to  have  general  depth 
profiles. 

Earlier  models  employed  dimensionless  ratios  for  sediment  properties,  and  sim¬ 
ilar  normalization  may  prove  convenient  in  applications  of  the  more  general  model. 
For  example,  sediment  mass  density  and  sound  speed  can  be  normalized  by  divid¬ 
ing  by  the  water  mass  density  and  water  sound  speed,  cw,  respectively.  Although 
density  normalization  is  mostly  a  matter  of  convenience,  Hamilton  and  Bachman  [7] 
argue  that  the  sound  speed  ratio  is  an  intrinsic  property,  insensitive  to  changes  in 
the  acoustic  properties  of  the  overlying  water.  Thus,  as  the  sound  speed  of  the  water 
at  a  given  site  undergoes  changes  due  to  warming,  cooling,  or  changes  in  salinity, 
these  ratios  remain  approximately  fixed.  Therefore,  one  need  change  only  the  water 
sound  speed  and  the  profile  for  sediment  sound  speed  will  automatically  adjust  to 
the  correct  value.  The  dimensionless  parameters  already  used  for  absorption  loss  8 
and  volume  scattering  <72  represent  an  attempt  to  remove  the  frequency  dependence 
of  these  parameters.  If  the  user  changes  the  frequency  but  leaves  all  other  param¬ 
eters  unchanged,  the  resulting  frequency  dependence  of  the  model  is  the  default 
dependence  that  we  expect  in  most  situations.  Of  course,  the  user  is  always  free  to 
specify  alternate  values,  including  a  nondefault  frequency  dependence. 

The  two  parameters  describing  bottom  roughness  are  defined  in  Section  3  and 
are  identical  to  those  employed  in  past  models  [2,  3].  In  previous  versions  of  this 
model,  the  default  value  7  =  3.25  was  recommended  when  a  measured  value  is  not 
available.  This  default  value  is  an  average  obtained  from  centimeter-scale  bottom 
roughness  spectra  obtained  by  stereophotography  [19]. 
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5  CONCLUSIONS  AND  RECOMMENDATIONS 

By  generalizing  the  propagation  calculations  within  a  prior  model  for  sediment- 
volume  and  interface  backscattering,  the  range  of  applicability  of  the  model  has  been 
greatly  increased.  More  general  profiles  for  sediment  acoustic  properties  now  per¬ 
mit  the  inclusion  of  small-scale  effects  such  as  loosening  of  the  upper  sediment  by 
bioturbation  and  such  large-scale  effects  as  consolidation  with  depth.  Although  the 
motivation  for  the  present  model  was  our  desire  to  join  earlier  models  covering  dis¬ 
joint  frequency  regimes,  we  went  beyond  merely  interpolating  between  the  two  prior 
models.  Our  generalized  treatment  of  sediment  propagation  has  introduced  new 
physics,  which,  in  some  cases,  leads  to  qualitatively  different  results.  Specifically, 
we  recognize  an  enhancement  in  scattering  strength  relative  to  the  gradient-free  case 
caused  by  the  presence  of  both  upward-  and  downward-propagating  sound  in  the 
sediment. 

While  the  model  is  capable  of  mimicking  the  angular  and  frequency  dependence 
seen  in  existing  backscattering  data,  it  has  not  been  validated  owing  to  a  lack  of 
measured  physical  parameters.  There  is  also  a  lack  of  wideband  backscattering  data 
spanning  low  to  high  frequencies  (hundreds  of  hertz  to  several  kilohertz).  In  the 
absence  of  such  data,  a  first  step  in  bringing  the  model  to  a  practically  useful  state 
is  to  constrain  the  sediment  profiles  based  on  current  geological  knowledge  and  to 
use  the  remaining  freedom  to  make  empirical  fits  to  existing  backscattering  data. 

Although  the  present  model  is  well  defined  mathematically,  and  we  have  eval¬ 
uated  it  for  some  nontrivial  cases,  care  was  continually  required  to  ensure  that 
particular  numerical  evaluations  always  yielded  correct  results.  Other  interested 
workers  will  have  to  exercise  similar  care  in  their  own  evaluations  of  this  model.  At 
present  we  could  not  supply  a  user-friendly  computer  implementation  that  would 
relieve  users  of  this  responsibility.  For  our  examples,  we  found  it  convenient  to 
choose  only  continuous  algebraic  profiles  of  sediment  properties.  Numerical  imple¬ 
mentations  working  from  tabular  profiles  and  allowing  for  discontinuously  layered 
sediments  would  clearly  be  desirable  as  well. 
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